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ABSTRACT 

We study some quantum mechanical aspects of dynamical black holes where 
the Vaidya metric is used as a model representing evaporating black holes. It 
is shown that in this model the Wheeler-DeWitt equation is solvable in whole 
region of spacetime, provided that one considers the ingoing (or outgoing) Vaidya 
metric and selects a suitable coordinate frame. This wave function has curious 
features in that near the curvature singularity it oscillates violently owing to large 
quantum effects while in the other regions the wave function exhibits a rather 
benign and completely regular behavior. The general formula concerning the black 
hole radiation, which reduces to the Hawking's semiclassical result when r = 2M is 
chosen, is derived by means of purely quantum mechanical approach. The present 
formulation can be applied essentially to any system with a spherically symmetric 
black hole in an arbitrary spacetime dimension. 
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Hawking's discovery that quantum black holes can evaporate by emitting ther- 
mal radiation [1] has stimulated not only active researches on quantum field theory 
in curved spacetime and a theory of quantum gravity but also a great deal of specu- 
lation about a grand synthesis of general relativity, thermodynamics, and quantum 
theory [2]. Unfortunately, however, in spite of much impressive effort it might be 
fair to say that in four dimensions little progress has been made in generalizing 
Hawking's "semiclassical" analysis to mathematically self-consistent, quantum me- 
chanically correct approach. 

Recently, we have developed a new approach which is purely quantum mechan- 
ical in that not only the matter fields but also the gravitational field are treated 
as quantum fields, in order to understand the black hole radiation [3, 4], and sub- 
sequently applied fruitfully to various problems relevant to quantum black holes, 
for example, the mass inflation in the Reissner-Nordstrom charged black hole [5], 
the three dimensional de Sitter black hole [6], the weak cosmic censorship [7], and 
the two dimensional dilaton gravity [8]. The key observation in this approach is 
that we can carry out the canonical quantization for a spherically symmetric sys- 
tem with the Vaidya metric representing the dynamical black holes if we limit our 
consideration to only the region near the apparent horizon and at the same time 
choose an appropriate coordinate frame. 

However, this approach gives us some frustrations. This is because within 
the approach it is impossible to gain useful informations about quantum struc- 
ture of spacetime in whole region through the wave function and it is also quite 
unclear what approximation method considering the specific region of spacetime 
corresponds to. 

The main motivation in this paper is to remove these disadvantages and put 
our approach on a more sound foundation. As a bonus, we will see that the 
wave function obtained in this way shows very interesting, physically reasonable 
features such that it fluctuates violently near the curvature singularity while it 
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behaves mildly in the other regions. In addition, we can derive the more general 
formula with respect to the mass loss rate of an evaporating black hole compared 
to the semiclassical approach. 

The action which we consider in this paper is of form 



where $ is a real scalar field. To show the four dimensional character explicitly 
we put the suffix (4) in front of the metric tensor and the curvature scalar. We 
follow the conventions adopted in the MTW textbook [9] and use the natural units 
G — h — c = 1. The Greek indices fj,,v,... take the values 0, 1, 2, and 3, while 
the Latin indices a,b, ... take the values and 1. Of course, the inclusion of other 
matter fields, the cosmological constants and the surface terms in this formalism 
is straightforward even if we confine ourselves to the action (1) for the sake of 
simplicity. 

After a general spherically symmetric reduction [10] 



S = I d^x^i% ( 



1 



(i) 
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ds 2 = ^g^dx^dx", 

= g ab {x c )dx a dx b + <p(x c f(d6 2 + sin 2 9d<p 2 ), 



(2) 



the action (1) can be rewritten as 




(3) 



Moreover, introducing the ADM parametrization 
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and n a which is unit vector normal to the foliations x° = const 



n a = (~, - — ), (5) 

a 07 



the action (3) becomes 

S = Jd 2 x[ X -a^ {1 - (n a d a cP) 2 + i(0') 2 - Kn a d a {4> 2 ) 
+ -<W 2 )} + 2 {(n^ a $) 2 - V) 2 } ] , 



(6) 



where K is the trace of the extrinsic curvature defined as y=<9a(v^~# n °) anc ^ 
^ = do and ^ = 9i are also denoted by an overdot and a prime, respectively. 

From (6), the canonical conjugate momenta can be read off 

P<t> = ~\fl n a d a (p - ^ K(f), 

Then the Hamiltonian can be found to be a linear combination of the Hamiltonian 
and the momentum constraints as follows: 

H = J dx l (aHo + pHi), (8) 
where the constraints are explicitly given by 



^ fifJ\2 2 V^ , 2 lVl 2 



(9) 
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#i = - + - 2p ; 7 - -p 7 y. (io) 

The canonical formalism explained so far gives a basis to construct a quantum 



theory of a system with a spherically symmetric black hole in the below. As a 
first crucial step toward the canonical quantization, we shall introduce the two 
dimensional coordinate x a by 1 



x a = (x°, x 1 ) — (v - r, r), 



(11) 



where v is the advanced time coordinate. Next let us fix the gauge symmetries by 
the gauge conditions 



where the black hole mass M is in general the function of the two dimensional 
coordinate x a . From (11) and (12) the two dimensional line element takes a form 



In previous works [3-8], the model which we have described till now and its 
variants have been widely used to perform the canonical quantization of dynamical 
black holes near the apparent horizon. The main new achievement in this article 
is to accomplish the quantization not only near the apparent horizon but over 
the whole spacetime region. This will enable the wave function to be explored in 

1 This choice of the coordinate system is very critical to reach the desired result. For in- 
stance, instead of (11) if we choose x a — (v,r) which reduces to the Eddington-Finkclstcin 
coordinate in vacuum, later analysis would reveal that the wave function is independent of 
the matter field and static. 




(12) 



ds 2 = g a bdx a dx b , 

= -(1 )dv 2 + 2dvdr. 



r 



(13) 
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considerable detail. To do so, since we are interested in the ingoing Vaidya metric 
[11], we impose the reasonable assumption on the dynamical fields 

$ = $(u), M = M(v), <p = r, (14) 

which is obviously consistent with the fields equations stemming from (1). Under 
the assumption (14), the canonical conjugate momenta (7) and the constraints (9) 
and (10) reduce to 



p$ = 2 d v $, 



1. 2M : 



2 



p <p =-d v M+- T -, (15) 



M 

7 



v^/Hq = 7 -Hi, 

1 2 2M 2 (16) 

where 7 = 1 + f rom (12). Here some comments are in order. First of all, the 
reason why the two constraints are proportional to each other in (16) is that there 
remains only one residual symmetry which generates the translation v — > v + e 
owing to the assumption (14) whose situation bears some resemblance to the case 
of the lightcone string field theory [12]. This corrects the wrong reasoning stated in 
previous our work [6] . Secondly, the dynamical degrees of freedom 7 corresponding 
to "graviton" are frozen because p 7 in (15) does not have the term including a 
differentiation with respect to v. This is nothing but a manifestation of the Birkoff 's 
theorem [9]. Finally, in (16), appears in the linear form so that this equation 
has the same structure as the Schrodinger equation if we regard a suitable function 
of (f> as the time in the superspace. 
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Substituting p$ = — i-g^ and = into the constraint (16) yields the 

Wheeler-DeWitt equation. 

. 1 d 2 d 2M 2 W 

(-ifim +l ~<M + — »*=°- (17) 

It is straightforward to obtain a special solution of (17) by the method of separation 
of variables. The result is 

$ = (£ e >/Z*(t;) + Ce -Vl*(«)) el ^#log 7) (18) 

where A, B, and C are integration constants. Let us examine the physical implica- 
tions of this wave function in detail. At the spatial infinity r — > oo, 7 — > 1 so that 
^ consists of only the ingoing matter field in the asymptotically flat spacetime as 
expected. And at the apparent horizon r = 2M(v), 7 = 2 so that \1/ is composed 
of the part of the ingoing matter field plus the part affected by the gravitational 
field there. Finally in the vicinity of the curvature singularity r — > 0, 7 — > 00 
thus \I/ oscillates violently. It is important to remember that the field 7 is exactly 
the dynamical degrees of the freedom of "graviton". Although classically these 
degrees of the freedom are killed owing to a spherically symmetric ansatz, it is of 
interest that they reappear and play an important role in quantum theory through 
the wave function. In other words, we can regard that strong quantum effects 
associated with the gravitational degrees of freedom cause the wave function to be 
singular at the singularity. Note that \I/ is completely regular over the whole region 
of spacetime except at the curvature singularity. Also notice that in deriving the 
Wheeler-DeWitt equation (17) and its solution (18) we have never appealed to any 
approximation method, thus they are essentially in the nonperturbative regime. 

Now under a general definition of the expectation value, the change rate of the 
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mass function of a black hole can be evaluated by means of (16) and (18) 

<8 V M>=-^. (19) 

If we set r to be equal to the black hole radius 2M , (19) reproduces the Hawking's 
semiclassical result up to a numerical constant [1]. But (19) is derived from the 
purely quantum mechanical arguments and provides the more general result than 
the semiclassical one although I have no idea how to determine the constant A 
within the present framework. In order to fix the numerical coefficient in this 
formula, maybe we would need to impose the boundary conditions on the wave 
function. Furthermore, this formula strongly suggests that the particles emitted 
by an evaporating black hole can be traced all the way back to not the surface 
of the horizon but the singularity. Incidentally, when we consider the outgoing 
Vaidya metric 

ds 2 = _(1 _ ^M)rf M 2 _ 2dudr (2Q) 

r 

instead of the ingoing one (13), we can repeat the previous arguments and it then 
turns out that we can reach similar results to the case of the ingoing Vaidya metric 
except a minor change where the coordinate v must be replaced by u. One weakness 
in the present model is the lack of outgoing radiation escaping to the future null 
infinity. It is valuable to mention that the loss of the black hole mass in (19) is 
controlled by sending in negative energy flux from the past null infinity, rather 
than by having positive energy flux radiated to the future null infinity. In future, 
it would be interesting to construct a model which takes account of both the energy 
fluxes at the same time. 

It is of interest to inquire whether the wave function has also a singular behavior 
near the the curvature singularity r = when we use a different kind of hypersur- 
faces, for example, x 1 = r = const which would be more appropriate inside the 
apparent horizon. Using the canonical formalism constructed previously [4] and 
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similar thoughts to the above, it is straightforward to derive the Wheeler-DeWitt 
equation 

l-?w -*>&-> v-tX*' ' (21) 

where this time 7 is given by 7 = —(1 — ^y-) > inside the apparent horizon. The 
wave function satisfying (21) is 

V]> = ^_g e -/A*(u) _|_ Q e -y/A${v)j e i(r-Mlogr)~^log 7 ^ ( 2 2) 

It is obvious that near the singularity the wave function (22) oscillates violently 
owing to the factors logr as well as 7 which would be again rooted to strong 
quantum effects of the gravitational degrees of freedom. It is valuable to comment 
two points here. One is that we have previously obtained the wave function with 
similar features to this in a different context where the radial quantization was 
carried out for "static" black hole models without matter field [13]. Their con- 
formity in the behavior of the wave function suggests that the singular behavior 
of the wave function at the singularity has a connection with quantum effects of 
the gravitational field. Of course, the quantum effects associated with the matter 
field would make the singular behavior there stronger. This point could be also 
understood from the fact that the change rate of the mass function (19) diverges 
at the singularity. The other is that we cannot connect smoothly the two wave 
functions (18) and (22) just on the apparent horizon. This is simply because we 
have used the different formulations. 

In conclusion, we have investigated some quantum aspects of a dynamical black 
hole corresponding to the Schwarzschild geometry in four dimensions. However, 
the present analysis is easily generalized to all spherically symmetric black holes 
in any spacetime dimension. One intriguing question in future is how the present 
analysis takes over if we use the Kruskal coordinate which enables us to slice the 
spacetime by the hypersurfaces extending from the left spatial infinity to the right 
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spatial infinity across the horizon and to find the wave function holding over the 
whole Kruskal extention of the Schwarzschild geometry. In this respect, Kuchar's 
work might be helpful [14]. 
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